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ABSTRACT 

A number of supernovae, classified as Type II, show remarkably peculiar properties 
such as an extremely low expansion velocity and an extraordinarily small amount of 
^^Ni in the ejecta. We present a joint analysis of the available observations for two of 
these peculiar Type II supernovae, SN 1997D and SN 1999br, using a comprehensive 
semi-analytic method that can reproduce the light curve and the evolution of the line 
velocity and continuum temperature. We find that these events are under-energetic 
with respect to a typical Type II supernova and that the inferred mass of the ejecta 
is relatively large. We discuss the possibility that these supernovae originate from the 
explosion of a massive progenitor in which the rate of early infall of stellar material 
on the collapsed core is large. Events of this type could form a black hole remnant, 
giving rise to significant fallback and late-time accretion. 

Key words: supernovae: general - supernovae: individual: SN 1997D - supernovae: 
individual: SN 1999br - methods: analytical 
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I INTRODUCTION 

Recently, a number of supernovae have been discovered that, 
according to their spectral properties, are classified as Type 

II but that, at the same time, show remarkably peculiar 
properties. The first clearly identified example was the ex- 
ceptionally faint SN 19970 in NGC 1536 (Turatto et al. 
( [l998| ); Benetti et al. ( |200l| )). SN 1997D showed both a very 
faint radioactive tail in the light curve, indicating an ejected 
^^Ni mass of only a few x 1Q~'^Mq, and a low expansion ve- 
locity of ~ 1000 km s_\ Modeling the spectra and light 



curve Turatto et al. (1998) concluded that a low energy ex- 
plosion in a 26 Mq progenitor star could successfully fit the 
early observations. The scenario of the low energy explosion 
of a high mass progenitor is also consistent with the late 
time ( ~ 400 days) spectral and p hotom etric data (Benetti 
et al. ( |200l| )). Chugai & Utrobin ( ^000| ) have presented an 
alternative analysis in which the progenitor was a star at 



Based on observations collected at the European Southern Ob- 
servatory and Cerro Totolo Inter- American Observatory, Chile 
(Proposals ESO 63.H-0141 and 66.D-0683). 



the low end of the mass range of core-collapse supernovae 
(8-12 Mq). 

The very low luminosity at discovery (10 times less than 
the peak luminosity reached by SN 1987A) and very small 
expansion velocity of the ejecta (3-4 times less than that of 
a normal Type II) of SN 1997D suggest that this explosion 
event was under-energetic. The mechanism that causes the 
energy of the explosion to be so low is a challenging and im- 
portant problem in the physics of core-collapse supernovae. 
In fact, the exact evolution of the star after neutrino reheat- 
ing depends on the rate of early infall of stellar material on 
the collapsed core and on the binding energy of the enve- 
lope. If both are large, as in high mass stars (M > 2OM0), 
the energy available to accelerate and heat up the ejecta can 
be greatly reduced and, eventually, may become insufficient 
to cause a successful explosion. Therefore, after the passage 
of the shock wave (and possibly the reverse shock formed at 
the H-He interface) a variable amount of matter may remain 
gravitationally bound to the c ollap sed remnant and fall back 
onto it (Woosley & Weaver (1995)). For this reason it was 
early suggested that SN 1997D could host a blac k hole rem- 
nant formed during the explosion (Turatto et al. (1998)) and 
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that the luminosity powered by fallback of envelope material 
onto the central black hole could emerge at about 10 00-12 00 
days after the explosion (Zampieri, Shapiro & Colpi (199S)). 
Unfortunately, to date it has not been possible to confirm 
or disprove observationally this prediction. 

After SN 1997D, a number of supernovae with similar 
observational properties have been identified. These objects 
appear to define a fairly hom ogeneous group of explosion 
events (Pastorello et al. ([2002D). As for SN 1997D, these su- 



pernovae provide a unique opportunity to probe the physics 
of the explosion and reach a better understanding of both 
the explosion mechanism and the conditions for the forma- 
tion of black hole remnants. In this Paper we present the 
results of a joint analysis of the observations of two peculiar 
Type II supernovae with very low luminosity and expansion 
velocity, SN 1997D and SN 1999br, which are representa- 
tive of the properties of the whole group. In particular, SN 
1999br represents the most extreme case of a low-luminosity 
event to date and appears to follow a recently reported 
correlation between expansion velocities of the ejecta and 
bolometri c lum inosities during the plateau phase (Hamuy 
& Pinto ( |2002| )V 

In Section H we describe the basic spectral and photo- 
metric data of these supernovae. Section ^ summarizes the 
semi-analytic method employed to model the light curve and 
the evolution of the line velocity and continuum tempera- 
ture of Type II supernovae. Finally, in section ^ the main 
results are presented and their consequences regarding the 
nature of the progenitors and the energy of the explosion 
are briefly discussed. 



2 PECULIAR TYPE II SUPERNOVAE 

Supernova 1997D, serendipitously discovered on 14 January 
1997 during an obser vatio n of the parent galaxy NGC 1536 
(De Mello & Benetti (1997)), is the first, clearly identified ex- 
ample of a peculiar Type II supernova with a very low lumi- 
nosity and expansion velocity (Turatto et al. (1995)). It was 
detected when it was already decaying from the plateau and 
was at least 2 mag fainter than a typical Type II supernova 
(Patat et al. (1994)). The decline rate of the last segment of 
the bolometric light curve is consistent with complete ther- 
malization of the gamma rays from the radioactive decay of 
''''Co into ^"^Fe. Assuming for SN 1997D the same deposition 
as in SN 1987A the ejected ^'^Ni mass is ~ 10"^ - lO^^M© 
(Benetti et al. (2001)). The spectra are dominated by a red 
continuum and P-Cygni profiles of H I, Ba II, Ca II, Na I 
and Sc II (see Figure hj). The most striking property of these 
spectra is the very low expansion velocity inferred from the 
spectral lines. The minima of the absorption lines give an 
expansion vel ocity of about 1100-1200 km s"^ for H (see 
Benetti et al. (2001) for details). 

Recently, other supernovae with properties similar to 
those of SN 1997D have been identified. A comprehensive 
analysis of the observations of these peculiar supernovae, 
including the recently discovered SN 2 001dc, will be pre- 
sented elsewhere (Pastorello et al. ( 2002| )). Here we focus on 
one particularly representative object, S N 19 99br in NGC 
4900, d iscov ered on 12 April 1999 (King (|l999|)). Filippenko 




et al. (199£) pointed out the exceptionally low luminosity 
of this event. Later, Patat et al. (1999) recognized that SN 
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Figure 1. Spectra of SN 1997D and SN 1999br at a phase ~ 100 
days. The spectrum of SN 1999br was obtained on 1999 July 20 at 
the ESO-Danish 1.54m telescope equipped with DFOSC (4- grism 
4). The spectrum of SN 1997D was taken on 1997 January 15 
with the ESO 1.52m telescope (B&C+grating 27). 



1999br shows other properties similar to those of SN 1997D, 
such as spectra with narrow P-Cygni lines, prominent Ba II 
features and a low continuum temperature. This supernova 
had a long plateau lasting at least 110 days with a mean 
luminosity of only <5 lO**^ erg s^^. It was extensively moni- 
tored during the first ~ 110 days with the CTIO and ESO 
telescopes, yielding extremely goo d pho tometric and spec- 
troscopic coverages (Hamuy et al. ( |2002| )). The spectrum of 
SN 1999br taken at a phase of ~ 100 days shows narrow 
metal lines (Figure 0). As in the case of SN 1997D, the red 
continuum and the strength of Ba II lines are caused by the 
low temperature of the ejecta. The U-B bands are strongly 
affected by line blanketing as the temperature decreases. 

As shown in Figure hi the spectrum of SN 1997D ob- 
tained soon after discovery and that of SN 1999br at ~ 100 
days are strikingly similar in both the continuum and the 
line components. Taking the ratio of the two spectra results 
in a rather flat curve where the main deviations are caused 
by slightly different line width (i.e. different expansion ve- 
locities of the ejecta) and spectral resolution. The expan- 
sion velocity of SN 1997D is only ~ 5% larger than that 
of SN 1999br. The close resemblance of these two spectra, 
the low inferred expansion velocity (~1000 km s^^) and the 
fact that the two supernovae have comparable luminosities 
strongly suggest that SN 1997D and SN 1999br may be sim- 
ilar events. In the following we will investigate in detail the 
consequences of this assumption. 

SN 1997D was discovered at the end of the plateau stage 
when the light curve was plummeting. The duration and lu- 
minosity of the plateau were inferred by comparison between 
models and observations and the consequent estimate of the 
explosion epoch ('-^ 6 days before discovery) was uncertain 
(Turatto et al. (1998)). On the other hand, the date of the 



explosion of SN 1999br has an uncertainty of only a few days 
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(thanks to a stringent pre-discovery limit). Then, from the 
similarity of the observational properties (luminosities and 
spectra) of SN 1997D and SN 1999br, we tentatively as- 
sume that the phase of SN 1997D at discovery is ~ 90 — 100 
days. This assumption will be checked by means of a detailed 
comparison of spectral-synthesis models with observations, 
presently under way. 



3 THE MODEL 

In order to determine the physical properties of a supernova, 
the observed light curve and spectra are compared to numer- 
ical simulations. Usually, it is important to carry out a pre- 
liminary investigation in order to obtain an approximate but 
reliable estimate of the physical conditions of the ejected gas 
and to establish a framework in which more detailed follow- 
up calculations can be performed. For these reasons, we have 
implemented a semi-analytic model in spherical symmetry 
that can provide a robust estimate of the parameters of the 
ejecta of Type 11 supernovae. The novelty of the present ap- 
proach lies mainly in the fact that the physical properties 
of the envelope are derived by performing a simultaneous 
comparison of the observed and simulated light curve, evo- 
lution of the line velocity and continuum temperature. The 
model follows the approach originally introd uced b y Arnett 
([1980|) and later developed by Arnett & Fu (|l989|) and Ar- 



nett (199C). The treatment of the motion of t he rec ombi- 
nation front follows in part the work of Popov (1992). The 



present analysis includes all the relevant energy sources pow- 
ering the supernova and evolves the envelope in 3 distinct 
phases that cover the whole evolution from the photospheric 
up to the late nebular stages. We include also the energy in- 
put from recombination that may be particularly important 
for low energy events like SN 1997D. The model has been 
tested against numerical radiation-hydrodynamic computa- 
tions in spherical symmetry under the same assumptions, 
giving good agreement. The details of the method and its 
valida tion w ill be presented in a companion paper (Zampieri 
et al. ( |2002| )). Here we briefly summarize the basic ideas. 

Full hydrodynamical calculations using realistic pre- 
supernova models show that the dynamical evolution of the 
envelope during shock passage is quite complex. The prop- 
agation of the shock determines how the explosion energy 
is distributed in the envelope of the progenitor star. The 
innermost layers comprised of helium and heavier elements 
transfer almost all of their kinetic energy and momentum 
to the Hydrogen envelope. The star mixes but does not ho- 
mogenize. The actual velocity, density and heavy element 
distributions of the post-shock material affect the light curve 
and estimation of the envelope mass in a major way. In the 
present analysis we do not consider the evolution of the star 
during this complex phase but rather assume idealized ini- 
tial conditions that provide an approximate description of 
the ejected material after shock (and possible reverse shock) 
passage, as derived from hydrodynamical calculations. The 
evolution starts at time tin (since core collapse) when the 
envelope is essentially free-coasting and in homologous ex- 
pansion. In realistic explosion calculations, it takes several 
hours in order for the envelope to relax to this state. We 
adopt tin — 5 hrs. At this stage, the velocity V of each gas 
shell is approximately constant and proportional to its po- 



sition r, V{r) — Vo{r/Ro), where Ro and Vq are the initial 
radius and velocity of the outermost shell of the envelope. 
Therefore, the radius of each gas shell increases linearly with 
time: for the outermost shell, R — Ro + Vo{t — tin)- The 
post-shock, ejected envelope is assumed to have spatially 
constant density p and total mass Menv = 4-KpR^/3. Mass 
conservation gives p = po{Ro/R)'^, where po is the initial 
density. In reality, the outer part of the star develops a steep 
power-law density structure that affects the light curve dur- 
ing the first few days after shock breakout (which is not 
included). However, this is only ~ 1% by mass of the star, 
while most of the stellar material resides in the inner part 
with roughly constant density. The latter region dominates 
the light curve after 10-20 days. The expansion velocity of 
the envelope as a function of interior (Lagrangian) mass 
m(r) = 47rpr^/3 is then V{m) = Vo(m/Me„„)^/^ (see Fig- 
ure B). The initial thermal-|-kinetic energy of the envelope 
is _E = (3/10) Menu Vo^//o, where /o is the initial fraction 
of kinetic energy. Elements are assumed to be completely 
mixed throughout the envelope. In our simple spherically 
symmetric model, their distribution depends only on r (or 
m). In particular. Hydrogen, Helium and Oxygen are as- 
sumed to be uniformly distributed, whereas ^^Ni is more 
centrally peaked (see Figure ^). Our assumptions about the 
velocity and elemental distributions provide only an approx- 
imate description of the actual post-shock structure of the 
envelope and should be regarded as a potential source of 
systematic uncertainty in the present model. 

The evolution of the supernova envelope starts at time 
tin and is schematically divided into 3 phases. During the 
first stage (a few tens of days) the envelope is hot and com- 
pletely ionized because of the energy deposited by the shock 
wave. In the following phase (from 30-40 days up to ~ 100 
days) , because of the decrease in temperature caused by ex- 
pansion and radiative diffusion, the envelope recombines and 
can be schematically divided into two regions, below and 
above the position of the recombination wavefront. In the 
third stage (after ~ 100 days) , the ejecta are completely re- 
combined and transparent to optical photons. The evolution 
is computed solving the energy balance equation for the en- 
velope gas. The thermal balance is governed by the competi- 
tion among the energy input from trapped gamma-rays, the 
PdV work and the energy losses through radiative diffusion. 
Assuming that radiation is in LTE with the gas throughout 
the envelope, the energy balance equation becomes a second 
order parti al differential equation for the temperature T. Ar- 
nett (1996) has shown that, under specific assumptions on 
the spatial distribution of ^^Ni, the energy equation (with 
appropriate boundary conditions) can be solved by sepa- 
ration of variables. The function that describes the radial 
dependence ip{r) is solution o f an e igenvalue equation (see 
e.g. equation [D.48] of Arnett (1996)). The final solution has 
the form 



V'(x)0(t), 



(1) 



where x — r/R, To is a constant reference temperature, 
(j}{t) is the function that describes the time dependence 
[4>{0) = 1] and, in the limit of zero mean free path, tp{x) = 
sin(7ra;)/(7ra;) (Arnett's "radiative zero" solution). The term 
(Ro/R)^ accounts for adiabatic 

radial temperature profile is T/To = 'tp{xy^* (see Figure B). 
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Figure 2. Initial (t = ti„) conditions of the post-shock envelope 
in our model. Upper panel: velocity V, density p and tempera- 
ture T distributions as a function of interior (Lagrangian) mass 
m. Bottom panel: Dependence of Hydrogen, Helium, Oxygen and 
Nickel abundances on m. 



Once a solution for is calculated, the emitted luminosity 
is approximately given by the expression for the diffusion lu- 
minosity L — — (47rr^c)/(3«:p)[9ar*/9r],.=jj, where k is the 
gas opacity. 

When the gas starts to recombine ai t = trec,Q, the 
part of the envelope below the position r; of the wavefront 
is assumed to be in LTE with radiation. We adopt the ap- 
proximation that in this region radiative diffusion can ef- 
fectively readjust the radial temperature distribution to the 
changing position of the outer b ound ary at ri = ri{t) ("slow 
approximation", Arnett & Fu (1989)). In this assumption, 
the spatial dependence of T is essentially given by Arnett 's 
"radiative zero" solution. Therefore, we search for approx- 
imate solutions of the energy equation of the form given 
in equation (|l|), with x replaced by x/ xj, Xj — Vi/R and 



ip{x/xi) — sin{nx/xi)/{ivx/xi) (Arnett (1996)). The func- 
tion that gives the time dependence of the temperature pro- 
file, 0(t), can be approximately computed setting the dif- 
fusion luminosity at recombination equal to the luminos- 
ity emitted by a blackbody at the effective temperature, 



(Popov ( |1992| )). This gives 



eff 
To 



,<hp_ 
dy 



(I) 



(2) 



where y = x/xi and T^n is the effective temperature, ap- 
proximately constant and nearly equal to the gas ionization 



temperature. Using equations (|l|) (with x replaced by x/xi) 
and (^), and integrating the energy equation over the ionized 
region, we obtain an equation for the motion of the recombi- 
nation wavefront Xi{t). Once a solution for Xi is computed, 
the bolometric luminosity from the inn er env elope is approx- 
imately given by (see Zampieri et al. (2002) for details) 



Lr,^L + A'Kr'iv,{aT^if/2 + pQ.on) , (3) 

where L = —{'iTvcaTQ)/{3K,po)[y'^dilj/dy]y^iRoXi(f>{t) is the 
diffusion luminosity, the second term on the right hand side 
is the total advection luminosity released because of the 
wavefront motion, Vi = XiR is the wavefront velocity relative 
to the envelope gas and Qion is the recombination energy per 
unit mass. 

In the recombined region {ri < r < R), the deposi- 
tion of gamma ray photon energy through Comptonization 
and photoelectric absorption of heavy elements is the dom- 
inant thermal and radiative process. The optical luminos- 
ity emitted in this region originates by the reprocessing of 
the gamma rays and becomes important only at the end of 
the recombination stage (that coincides with the end of the 
plateau), when the radioactive decay time of ^^Co becomes 
comparable to the expansion timescale. An approximate ex- 
pression for this luminosity can be obtained neglecting the 
internal energy and the PdV work in the energy equation 
and integrating it over the interval ri < r < R. The total 
bolometric luminosity of the envelope during this phase is 
then 



(4) 



Lr, + 4-KpoRoXN^,ofit) / X^Ni{x)dx 



where Xjvi,o and (,Ni{x) are the central mass fraction and 
radial distribution of ''^Ni, f{t) = [3.9 x loi°e-*/^"' + 
7.2 X 10^(e-*''^co _ e"*/^"')] erg g~^ s-\ rm = 8.8 days 
and TCo = 111 days a re th e Nickel and Cobalt decay 
times. Following Arnett ( 1996 ), we take ^Ni{x) = tp{x) = 
sin(7rx)/(7ra::) (see Figure g, where the ^®Ni abundance is 
shown as a function of m). 

In the third stage, when the envelope has completely 
recombined, the ^''Co radioactive decay time is shorter than 
the expansion timescale and the luminosity is given by the 
second term in equation (m with Xi = 0. 



4 RESULTS AND DISCUSSION 

Following the discussion in § 2, we assume that the phase 
of SN 1997D at discovery is ~ 90 - 100 days and that the 
early light curve and the evolution of the line velocity re- 
semble closely those of SN 1999br. Figure ^ shows the evo- 
lution of the UBVRI luminosity, the velocity of metal (Sc 
II) lines and the continuum temperature of SN 1997D and 
SN 1999br along with the results of our semi-analytic model 
calculations. In the top panel, the observed and calculated 
light curves of SN 1987A are also shown for comparison. 
The continuum temperatures have been estimated fitting 
the spectra with a Planckian and are affected by a significant 
uncertainty (±500 K) because of the severe line blanketing 
at short wavelengths (< 4500 A, see Figure |l|). The line ve- 
locity and continuum temperature inferred from the model 
refer to the photospheric epoch. Hence the last two plots 
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Figure 3. UBVRI luminosity, L, velocity of metal (Sc II) lines, 
Vg^i, and continuum temperature, Tc, as a function of time for 
SN 1997D (circles) and SN 1999br (triangles). The light curve 
of SN 1987A is also shown for comparison. The adopted dis- 
tance modulus and the estimated total reddening are /i = 18.49, 
E{B -V) = 0.18 for SN 1987A, = 31.29, E{B - V) = 0.02 
for SN 1997D and = 31.19, E{B - V) = 0.02 for SN 1999br 
(see Pastorello et al. (2002) for details). The solid, long-dashed 
and short-dashed lines represent the model curves. The veloci- 
ties in the mid panel are those of the gas at the position of the 
recombination wavefront Vg^i = XiVo. The asterisk is the contin- 
uum temperature inferred from the spectral synthesis model of 
SN 1997D at discovery (Turatto et al. (1998)) moved to a phase 
of ^ 100 days. 



in Figure ^ are truncated at the end of the recombination 
phase (~ 110 days). 

Considering the approximations adopted in the model, 
the general agreement with observations is satisfactory. We 
emphasize that obtaining a simultaneous "fit" of the light 
curve and the evolution of the line velocity and continuum 
temperature is an essential requirement to obtain a mean- 
ingful estimate of the parameters of the ejected envelope. 
In particular, the long plateau, the apparent break in the 
line velocity profile at ~ 30 days and its fast decline impose 
rather severe constraints on the model. The break appears 



to be related to the onset of recombination, that terminates 
when the light curves plummets at 110 days. 

The parameters of the post-shock, ejected envelope re- 
quired to reproduce the observations of SN 1997D and SN 
1999br are listed in Table |^. Only the values of the radius 
Ro, mass Men., "'Ni mass Mmz and outer velocity Vo (and 
hence initial energy) of the envelope have been significantly 
varied. The other model parameters were maintained essen- 
tially equal to the values required to reproduce the observa- 
tions of SN 1987A. In particular, a colour correction factor 
fc = Tc/Tfijj = 1.1 — 1.3 has been adopted to reproduce 
the observed continuum temperature for all the supernovae. 
Small colour corrections are rather common in stellar and 
supernova atmospheres and are induced by distortions of the 
continuum caused by radiative transfer processes (e.g. scat- 
tering, wavelength d epend ent opacities; see e.g. Eastman, 
Schmidt & Kirshner ( |l99^ )). 

Interpolating simultaneously the light curve and the 
evolution of the line velocity and continuum temperature 
results in a rather robust "fit". The estimated value of the 
initial thermal-|-kinetic energy of the ejecta E indicates that 
both events were rather under-energetic compared to a typi- 
cal Type II supernova. The inferred ^'^Ni mass of SN 1999br 
is extremely small, testifying that the energy available to 
produce and eject nucleosynthetic elements was very low. 
Furthermore, because of the low luminosity in the plateau 
stage, the post-explosion envelope is rather compact. The 
ejected envelope masses are quite large, comparable to those 
required to reproduce the plateau of typical Type II super- 
novae. In particular, the ejected envelope mass of SN 1997D 
is almost three times larger than that estimated by Chugai & 
Utrobin (2000) and only 30% smaller than that inferred by 
It is worth noting that the gross prop- 



1998 



Turatto et al. ( 

erties of the light curve and line velocity of SN 1997D and 
SN 1999br can be roughly accounted for by SN 1987A-like 
parameters, but simply decreasing the expansion velocity Vb 
(and hence the energy E) and ^^Ni mass. 

Although it is not straightforward to determine the er- 
ror in the parameters, we estimate that the intrinsic uncer- 
tainty ofthe "fit" is not larger than ~ 30% (see also Zampieri 
et al. (2002)). Additional sources of systematic errors are re- 
lated to the approximations introduced in the present anal- 
ysis, in particular to the choice of the initial conditions. Al- 
though it is known that signif icant mixing occurred in SN 
1987A (see e.g. Woosley (1988)), this effect may be less pro- 
nounced in other supernovae. The light curve is very sensi- 
tive to the prescription for mixing and to the actual velocity 
distribution as a function of mass V{m). In particular, the 
luminosity and duration of the plateau depend on the energy 
and mass of the high velocity. Hydrogen rich part of the en- 
velope. If the innermost Helium and heavier elements layers 
did not mix appreciably with the Hydrogen envelope, they 
would not produce any observable effect (having very low 
energy) . In this case, the estimated value of M^nv would sim- 
ply refer to the Hydrogen envelope mass. Then, in general, 
Menv represents a lower limit to the total meiss of the ejecta. 
Furthermore, the estimate of the envelope mass is sensitive 
to uncertainties in the actual value of the gas opacity and 
the details of the recombination physics. Here we adopted 
K = 0.2 cm^ g"\ trecfi ^ 30 days and T^// = 4000 - 4400 
K, values similar to those used for the "fit" of SN 1987A. 
Only a larger opacity and/or a delayed onset of recombina- 
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Table 1. Parameters of the semi-analytic model 





Ro 






Vb 


E 


/o 


K 


^rec.O 




/c 




(IQi^ cm) 


(Mo) 


(Mq) 


(10* cm s-i) 


(10^1 erg) 




(cm^ g^"*") 


(days) 


(K) 




SN 1987A 


5 


18 


7.5 X 10-2 


2.7 


1.6 


0.5 


0.2 


25 


4800 


1.1 


SN 1997D 


9 


17 


8 X 10~3 


2.1 


0.9 


0.5 


0.2 


30 


4400 


1.2 


SN 1999br 


7.5 


14 


2 X 10-3 


1.9 


0.6 


0.5 


0.2 


30 


4000 


1.3 



Ro is the initial radius of the envelope 
Menv is the ejected envelope mass 
Mjvi is the mass of ^^Ni 

Vb is the velocity of the envelope material at the outer shell 

E is the initial thermal -|-kinetic energy of the ejecta 

fo is the fraction of the initial energy that goes into kinetic energy 

K is the gas opacity 

trec,o is the time when the envelope starts to recombine 
Ti,f f is the effective temperature during recombination 
fc = T/T^ff is the colour correction factor 



tion could in principle allow for smaller envelope masses, but 
this is not in agreement with what found f or SN 1987A and 
other Type II supernovae (Zampieri et al. (2002)). The fact 
that the line velocities at ~ 100 days are slightly larger than 
what predicted by the model (see Figure ^ seems to indicate 
that the assumption of spatially constant density is only ap- 
proximately correct and that more mass is concentrated in 
the innermost, low velocity part of the envelope. This effect 
may reduce the estimated energy E because the bulk of the 
kinetic energy is carried by the outer, high velocity layers. 

We emphasize that the present analysis for SN 1997D 
is based on the assumption that it is an event similar to SN 
1999br. However, although the inferred parameters of the 
ejecta of SN 1997D rely on this hypothesis, the estimates for 
SN 1999br are certainly valid because they are not affected 
by uncertainties on the early light curve, the duration of the 
plateau and the evolution of the line velocity. 

Determining the mass of the progenitor is quite a diffi- 
cult task and, unless a pre-dis cover y identification is avail- 
able (see, e.g., Smartt et al. (2002)), the inferred value is 
usually rather uncertain. Adding to the mass of the envelope 
reported in Table ^ the mass of the collapsed core (« 2Mq), 
the progenitors of SN 1997D and SN 1999br have at least 19 
and 16 Mq , respectively. We stress that these are most prob- 
ably lower bounds and the actual values of the progenitor 
masses are likely to be larger. Therefore, the present esti- 
mates situate SN 1997D and SN 1999br in the intermediate 
mass progenitor range and rule then out an origin from the 
low end of the mass range of core-collapse supernovae. 

It should be noted that the results of the present anal- 
ysis are essentially insensitive to the inclusion of Rayleigh 
scattering from neutral Hydrogen. In fact, towards the end 
of the plateau, Tra ~ Krapr ^ 0.1, where r ~ 10^^ cm, 
p ~ IQ-^^ g cm-'^ and the Rayleigh scattering opacity 
Kra ~ IQ-* cm^ g-^. Therefore, the outer recombined re- 
gion is essentially transparent to Rayleigh scattering. 

Understanding why the energy of these peculiar Type II 
supernovae is low is of paramount importance in connection 
with the hydrodynamics of the explosion and the formation 
of the central compact object. Because the mass of the iron 
core before the explosion does not vary signifi cantly with 
progenitor mass Af, (Woosley & Weaver ( 1995 )), the gravi- 



tational potential energy liberated during the collapse of the 
core is roughly independent of M, . Therefore, the energy of 
the ejecta depends mostly on neutrino reheating and the 
hydrodynamic interaction of the supernova shock with the 
surrounding material. To unbind a spherically symmetric 
envelope, the shock must overcome the ram pressure of the 
gas that st arted to accrete after the collapse of the core (see, 
e.g.. Fryer (1999)). With increasing progenitor mass both the 



ram pressure and the binding energy of the envelope become 
larger. Therefore, comparatively less energy remains avail- 
able to heat up and accelerate the ejecta. The final luminos- 
ity and expansion velocity are then small and the resulting 
supernova is under-energetic (Zampieri ( |2002| )). In this re- 
spect, it is tempting to note that the difference between 
the inferred total (kinetic-|-thermal-|-binding) energy of SN 
1987A (~ 2.5 X lO'^^ erg) and the ejected kinetic-|-thermal 
energy of SN 1997D is ~ 1.5 x lO^'^ erg, comparable to the 
binding energy of the en velope in a 25 — 30Mq progenitor 
star (Woosley & Weaver ( |l995| )). Therefore, in SN 1997D, a 
large fraction of the gravitational potential energy liberated 
by the collapse of the core may have been spent in trying to 
unbind the massive envelope. 

An important consequence of this scenario is that, be- 
cause of the low kinetic energy acquired by the envelope, a 
large fraction of the stellar material (in particular the inner- 
most layers that have smaller velocities) is likely to remain 
gravitationally bound to the core after shock passage and 
fall back onto it. This is confirmed also by the small amount 
of ^^Co present in the ejecta. The fallback of stellar ma- 
terial may also turn the newly formed neutron star into a 
black hole. If this happens, the late time fallback onto the 
central black hole may give rise to detectable emission of 
radiati on wi th a characteristi c pow er-law decay (Zampieri 
et al. ( [1998| ), Balberg et al. ( ^OOOj )). The detection of this 
emission in the late time light curve would provide the first 
direct evidence for the presence of a black hole in the site of 
its formation. 

We note that the overall picture may be different for 
a supernova from a massive progenitor whose post-shock 
envelope retains a significant fraction of its initial angular 
momentum. In this case, the hydrodynamics of the explosion 
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becomes more complex and may give rise to a high e nergy, 
asymmetric explosion (MacFadyen & Woosley (1999)). 

The present analysis indicates that the parameters of 
the ejecta of SN 1997D and SN 1999br are consistent with 
them being intermediate mass core-collapse supernovae. If 
the mass of their progenitors is sufflciently large, they could 
form a black hole remnant, giving rise to significant fallback 
and late-time accretion onto the central compact object. In 
order to confirm these findings, more detailed investigations 
using radiation hydrodynamic simulations and spectral syn- 
thesis calculations are being planned. Clearly, monitoring 
this type of supernovae from discovery to late phases is of 
the utmost importance. 
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